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SUMMARY 


Shuttle flight test data have been used to determine values for the 
short-period parameters. The best identified, as judged by its estimated standard 
deviation, was the elevon effectiveness parameter C,,^. However, the scatter 

about the preflight prediction of Qngg was large. Other investigators have 
suggested that adding nonlinear terms to the mathematical model used to identify 
Cn^e could reduce the scatter. 


The results of this investigation show that is the only identifiable 

nonlinear parameter applicable and that the changes in Gmgg values when 


" m 6e2 


is included are in the order of ten percent for the data estimated. 


INTRODUCTION 

The longitudinal parameters that describe up to 90 percent of the short-period 
motion of the shuttle vehicle have been determined using flight test data. These 
parameters are shown plotted versus Mach number in figure 1 taken from reference 1. 
Of the parameters shown, C^g has the smallest standard deviations and the 

largest sensitivities, and it was considered the parameter best identified. The 
trends of the parameter with Mach number matched those of the preflight predictions; 

however, the scatter of the values extracted was large and reduced the confidence in 

the actual values determined. 

Investigators at Johnson Space Center have found that by using a nonlinear 
model and including the term Cn,^ e2 , the identif lability of Cm 6e was 

improved, and its value was close to that of the preflight predictions. The purpose 

of this paper is to use Shuttle Discovery flight test data to determine the 

identif lability of the Cm^ e 2 parameter, to check the possibility of including 

other nonlinear terms in the model, and to compare the values of Cmg e determined 
using the linear and nonlinear extraction models. 
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SYMBOLS 

acceleration measured along the Z body axis, g units 

wing mean geometric chord, m (ft) 
pitching-moment coefficient, M^/qSc 

axial-force coefficient, F^/qS 

vertical-force coefficient, F^/qS 

force along X and Z body axes, respectively, N (lb) 


acceleration due to gravity, m/sec2 (ft/sec2) 
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moment of inertia about X, Y, and Z body axes, 
respectively, kg-m2 (slug-ft2) 

product of inertia, kg-m2 (slug-ft2) 

pitching moments, N-m (ft-lb) 

mass, kg (slugs) 

rate of roll, pitch and yaw, rad/sec or degrees/sec 
dynamic pressure, N/m2 (slug/ft2) 
wing area, m2 (ft2) 

velocity along X, Y, and Z body axes, respectively, m/sec 
( ft/sec) 

airplane total velocity, m/sec (ft/sec) 

angle of attack, rad or degrees 

elevator, deflection, rad or degrees 

pitch angle, roll angle, rad or degrees 
air density, kg/m^ (slug/ft3) 


Stability and control derivatives referenced to a system of body axes with the 


origin 

i at the 

airplane center 

of gravity: 
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MODEL VERIFICATION 

In order to verify that the assumed nonlinear model best describes the vehicle 
motion several questions must be answered: 
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1. Is Cm, identifiable? 

“oe 

2. Does the addition of (^2 improve the identif lability of the other 
parameters in the model, especially C ni) j e ? 

3. Is C,^ 2 the best parameter for an extension of the mathematical model with 
linear aerodynamic parameters? 

4. Do the nonlinear terms added make sense physically? 

5. Does the nonlinear aerodynamic model have better prediction capabilities than 
the linear one? 

In the discussion that follows these questions will be answered. 


EQUATIONS 


The equations used in this study are standard equations used to describe the 
longitudinal motion of a vehicle in the atmosphere. For this study, the aerodynamic 
model equations are made nonlinear by the addition of terms such as 

( " m a6e ot ^ e or c n> a 2 a ^* The basic equations can be found in many references 

such as reference 2 and are repeated here for convenience. 
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u = -qw + rv - gsin0 + — [ C + C (a - a ) 1 

2 m L x- x t 1 

0 a 


w = -pv + qu + gcos9 cosiji + — f C + C_(a-a) + C_ (6 - 6 ) 1 

2 m 1 Z_ 7. t Z.e e^ J 
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DISCUSSION 

Longitudinal maneuvers from the first two Shuttle Discovery flights were 
examined using a maximum likelihood parameter extraction algorithm (ref. 3). 

Thirteen runs were selected covering a Mach number range from 22 to 1. The results 
of processing these data for several assumed mathematical models are shown as figure 
2. The values determined for g 2 and Cm^g at various Mach numbers are 

shown in Table I. The addition of the C m ^ e 2 term seemed to have the greatest 
effect on the values determined for the other parameters at the highest Mach 
numbers. The C,,^ parameter had the greatest variations when C m g e 2 was added 
to the mathematical model. The C,,,^ parameter also showed some change when 
C m g e 2 was added to the model, but below Mach 15, this change was generally 10 
percent or less. 
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The question of identif lability can be discussed by considering Table II. The 
values of sensitivity and standard deviation given are typical of the values seen 
for all thirteen runs and are a reasonable basis for the discussion of identif labil- 
ity. The sensitivity referred to is the variation in the output states with respect 
to a perturbation in a specified parameter assuming the other parameters are fixed. 
The greater the sensitivity the more identifiable the parameter. An examination of 
the sensitivities implies the C,n§ e 2 * s * ess identifiable than Cn,g e , as 

identifiable as Cn^ , and more identifiable than Cz a . The estimated standard 
deviations indicate that the values for were well determined since the 

estimated standard deviation was less than one-tenth of the extracted value. Based 
on this assessment, Qn$ e 2 was considered identifiable. In general, the addition 
of Cmg e 2 resulted in Cmgg variations of about 10 percent. 

The addition of C m ^ e 2 to the model appeared to improve the identif lability 

of C m ^ e . The sensitivities increased, and the estimated standard deviations as 

a percent of the extracted parameter did not degrade, implying an improved identif i- 
ability of C^g. Also, the change in the values determined for were 

generally in a direction so that they were closer to the preflight estimates, 
although this improvement was small for most runs. 

Since the shuttle maneuvers resulted in small amplitude responses, the values 
of even the best determined parameters show considerable scatter. In particular, in 
figure 1 was well determined for all runs, but showed significant run-to- 

run scatter. In this case, when a system is poorly excited, any additional param- 
eter can many times improve fit and identif iability (see ref. 4). To demonstrate 
that C m g e 2 was the best parameter to add to the mathematical model, several 

other parameters were tried in the model, specifically, Cm ct g e and Cm a 2* The 
effect on C m g e of adding Cp^g is shown in figure 2 for several runs. As 
can be seen, the value of Cp^g was not significantly changed when Qn^g was 

added to the model. Also, in the course of this analysis it was found that 
^ m a6e P rovec * to tre much less identifiable than C m g e 2 when the sensitivity 

and standard deviations were compared for the same flight data run with Cp,g e 2 

replaced in the model by Cp^g. The parameter Cpi a 2 proved to be totally 

unidentifiable and resulted in no changes in any of the parameter values. The 

improvement seen when C m< j e 2 was added to the model is not an artifact of the 

identification procedure but the legitimate result of adding a term that actually 
describes part of the vehicle motion, as evidenced by the identif iability discussed 
earlier in this section. 

The next question that must be answered is whether or not the parameter to be 
added makes sense physically. The values determined for Cmg e 2 were negative 

(see Table II). This implies a C m versus 6 e relation as shown in figure 3. The 

result of this relation would be that a negative 6e would result in less response 
than a positive 6 e . Figures 4 and 5 show that this is indeed the case for the 


4 



shuttle vehicle. For elevon inputs that are essentially equal in magnitude but 
opposite in sign, the figures show that the responses are at least 40 percent 
greater for the positive 6e input. 

Finally, since C n)( j e 2 is the only nonlinear term identifiable, is the 
predictive capability of the model including C ni( j e 7 better than the model where 
C m g e 2 is not included? Figures 6 and 7 show that using the nonlinear model 
including C m ^ e 2 does result in a slightly better fit to a set of responses that 
were not used to determine parameter values. 


CONCLUDING REMARKS 

The inclusion of 0,^2 in the mathematical model describing the 

longitudinal motions of the shuttle vehicle was evaluated using flight test data 
from the Discovery vehicle. The parameter was found to be identifiable, 

and the values determined were well defined. The parameters can be justified physi- 
cally by considering the increased responses of the vehicle to a positive elevon 
deflection when compared to the response from a negative deflection. Inclusion of 
the Cmg e 2 parameter caused C m ^ e values to move closer to the preflight 

predictions, but the changes seen were small. The model that included C m( j e 2 was 

used to predict the response of a run not used for identification, and the fit was 
slightly better with C ni( j e 2 in the model then when C m g e 2 was not included. 

These results imply that incorporating C nt g e 2 produces an improved model for the 

identification of the longitudinal aerodynamic parameters of the shuttle vehicle. 
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TABLE I 

C AND C VALDES AT 

m 6e 2 m 6e 

SPECIFIC MACH NUMBERS 


C „ 
m 6e 2 

C 

Mach 

m 6e 

22 

-.34 

-.170 

18 

-.35 

-.176 

18 

-.69 

-.235 

11 

-.26 

-.16 

8 

-.315 

-.142 

8 

-.32 

-.17 

6 

-.21 

-.117 

5.5 

-.21 

-.110 

4 

-.13 

-.087 

4 

-.09 

-.087 

1.7 

-.57 

-.15 

1.6 

-.051 

-.17 

1.0 

-1 .05 

-.42 
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TABLE II SENSITIVITIES AND STANDARD DEVIATIONS FOR SELECTED SHUTTLE RUNS 
WITH C m5e 2 AND WITHOUT C^ e2 IN THE EXTRACTION MODEL 


Mach = 8 


Flight Number = 19 


Parameter Value Standard Sensitivity Value Standard Sensitivity 
Deviation Deviation 


C Tm$e 

^m <$ e 2 





.0012 

.15 


. 19E+8 

.22E+8 


.004 

•85E+4 

.0008 

•52E+5 



Mach = 8 


Flight Number = 14 




.018 

. 1 5E+4 

.0038 

. 1 7E+5 



Standard 

Deviation 

Sensitivity 

.0058 

.2E+6 

.62 

.34E+3 

.0085 

.96E+4 

.0019 

.67E+5 

.016 

.8E+4 
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TABLE III SENSITIVITIES AND STANDARD DEVIATIONS FOR SELECTED SHUTTLE RUNS 

UTTD n Aim UTWIIAim O r«» »■■■■ ■ ww m* » \ 

"Am ^mg e 2 "Ainvui ^m,s e 2 xw AnL uiKAtnun nuutL ^u/niifiuku; 


Mac! 

■ ■■■ ... I. 1 

i = 6 


Flight Number 

= 19 

Parameter 

Value 

Standard 

Deviation 

Sensitivity 

Value 

Standard 

Deviation 

Sensitivity 


.13E+9 

.42E+4 


.00066 

.11 


. 15E+9 
.98E+4 


mo g 

C n>6 ( >2 


.0013 


.00035 


.85E+4 


. 16E+6 


.002 

.37 E+4 

.00036 

.011 

.37E+6 

.15E+6 


Mach = 6 


Flight Number = 14 


Parameter Value Standard Sensitivity Value Standard Sensitivity 
Deviation Deviation 


-.033 

1.28 

-.715 

-4.69 

-.066 

.0004 

-.20 


.0014 

.23 


•21E+7 

.28E+4 


•52E+4 


. 19E+7 
.25E+4 


. 4 1 E+4 


«>o e 
c m6 p 2 


-.095 .0025 


.33F.+5 


. 10E+6 
.84E+4 
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(a) Change In normal force with alpha 
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(b) Static stability 
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Figure 1. - Dominant longitudinal derivatives affecting the short period 
motion of the Shuttle versus Mach number (derivatives per 
degree) . 
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Figure 2. - Dominant longitudinal derivatives affecting the short period 

motion of the Shuttle versus Mach number from Discovery flight 
data and using three mathematical models. 
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Figure 3. 


versus 6 e curves for linear model and non-linear model 
including C mfie 8 e + Cn,^ fi e 2 terms 
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Time, sec 

Figure 4* - Time history of doublet maneuver with elevon initially negative 
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Figure 5. - Time history of doublet maneuver with elevon deflection 
initially positive. 
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Figure 6. - Prediction of vehicle motion using linear model (only Cm 
included in model). 
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